GENERALIZED FOSTER'S IDENTITIES 
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Abstract. Foster's network theorems and their extensions to higher orders involve re- 
sistance values and conductances. We establish identities concerning voltage values and 
conductances. Our identities are analogous to the extended Foster's identities. 



1. Introduction 

Foster's first identity, which is also known as Foster's network theorem, was initially proved 
by R. M. Foster [Folj in 1949 by using Kirchhoff 's rule. In 1990, D. Copersmith, P. Doyle, P. 
Raghavan, and M. Snir [CDRSJ gave another proof of this theorem by using random walks 
on graphs. Independently, in 1991, P. Tetali [TP1] gave a third proof of this theorem again 
by using random walks on graphs. (See also [Bj3[ Theorem 25 and exercise 23 in Chapter 
IX]). Later in 1994, Tetali |TP2j gave a fourth probabilistic proof by using an elementary 
identity for ergodic Markov chains. 

In the late 1980's, T. Chinburg and R. Rumely introduced a canonical measure fi can on 
a metrized graph. As a consequence of the fact that n can has total mass 1, Chinburg and 
Rumely [CR, remark on pg. 26] obtained an identity. In 2003, M. Baker and X. Faber [BF, 
Corollary 6] showed that this identity is equivalent to the Foster's first identity. This became 
another interesting proof of the Foster's first identity. 

In another direction, L. Weinberg [W] in 1958, D. J. Klein and M. Randic |D-Ml Corollary 
CI] in 1993, and R. Bapat jRBll Lemma 2] in 2004 proved identities that are equivalent to 
the Foster's first identity. The properties of discrete Laplacian and its pseudo inverse were 
used in the articles |D-M] and [RBI] . 

In 1961, Foster [Fo2] proved another identity which we call Foster's second identity. In 
2002, J. L. Palacios [P] gave a probabilistic proof of the Foster's second identity, generalizing 
the arguments used in |TP2j . He also showed how to extend the Foster's identities to graph 
paths with more edges and he proved a third identity. The Foster's first and second identities, 
and the third identity due to Palacios are about the sums of the resistance values over the 
graph paths consisting of one, two and three edges, respectively. 

In 2007, we [Cll Section 5.5] proved four identities that involve voltage values rather 
than resistance values. Our method was to use the properties of discrete laplacian and its 
pseudo inverse. As a corollary to these identities, we obtained Foster's 1st, 2nd, 3rd, and 
4th identities. 

In 2008, E. Bendito, A. Carmona, A. M. Encinas and J. M. Gesto jBCECl Proposition 
2.3] extended the Foster's identities to higher orders. Their formula contains all previously 
known Foster's identities as specific cases. 
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In this paper, we extend the identities concerning voltage values \C1\ Section 5.5] to higher 
orders. This can be found in Theorem 14. 21 which is our main result. In this way, we generalize 
the extended Foster's identities. We provide two proofs for Theorem 14.21 Both of the proofs 
rely on the properties of the discrete Laplacian and its pseudo inverse. In the first proof, 
we use a remarkable relation between the resistance values and the pseudo inverse L + of the 
discrete Laplacian [RBI J [RB2J (see also Lemma [3T4|) to express the voltage values by using 
L + (see Lemma [3 .51) . Similarly, in the second proof, we use an equally interesting relation 
between the resistance values and some "equilibrium measures" z/'s [BCEGl Proposition 2.1] 
to express the voltage values by using z/'s (see Proposition |4.4p . More information about the 
equilibrium measure can be found in the article |BCEG| and the related references therein. 

Note that there is a 1 — 1 correspondence between the metrized graphs and the equivalence 
class of finite connected weighted graphs in which the weight of an edge is the reciprocal of 
its length [BRl Lemma 2.2]. We can also view such a graph as a resistive electrical network 
in which the resistance of each edge is the same as its length. We will work with metrized 
graphs in this paper. However, the results and their proofs are valid for the corresponding 
weighted graphs and the resistive electrical networks. 



2. The voltage and resistance functions 

A metrized graph T is a finite connected graph equipped with a distinguished parametriza- 
tion of each of its edges. One can find other definitions of metrized graphs in the articles 
[CR], [BRj, [ZhJ Appendix] and [BF]. 

A metrized graph can have multiple edges and self- loops. For any given p G T, the number 
of directions emanating from p will be called the valence of p, and will be denoted by v(p). 
By definition, there can be only finitely many p G T with v(p) ^ 2. 

For a metrized graph T, we will denote its set of vertices by V(r). We require that V(T) 
be finite and non-empty and that p G V(T) for each p G T if v(jp) ^ 2. For a given metrized 
graph T, it is possible to enlarge the vertex set V(r) by considering more additional points 
of valence 2 as vertices. 

For a given metrized graph T, the set of edges of V is the set of closed line segments with 
end points in V(r). We will denote the set of edges of Y by E{Y). We will denote #(V(r)) 
and #(E(T)) by n and e, respectively if there is no danger of confusion. We denote the 
length of an edge G E(T) by Lj. 

In the article |CR| . the voltage function j z (x,y) is defined and studied as a function of 
x,y, z G T. For fixed z and y it has the following physical interpretation: when V is viewed as 
a resistive electric circuit with terminals at z and y, with the resistance in each edge given by 
its length, then j z (x,y) is the voltage difference between x and z, when unit current enters 
at y and exits at z (with reference voltage at z). 

For any x, y, z in T, the voltage function j x {y,z) on T is a symmetric function in y 
and z, i.e., j x (y,z) = j x (z,y). Moreover, for any x, y, z in T it satisfies j x (x,z) = and 
jx{y, y) = r(x, y), where r(x, y) is the resistance function on T (for more information see the 
articles [CR] and |BR[ sec 1.5 and sec 6], and [Zh, Appendix]). 

For any x G T, by circuit theory we can transform T to an V-shaped graph with the same 
resistances between x, p, and q as in F (see |C2[ Section 2] for more details). This is shown 
in Figure dj with the corresponding voltage values on each segment. Therefore, 
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Figure 1. Circuit reduction with reference to 3 points x, p and q. 

1 r(p,x) = j p (x,q) +j x (p,q), r(q,x) = j q (x,p) + j x (p,q), r(p,q) = j q (x,p) + j p (x,q). 
Then it follows from (TjQ) that 

(2) 2j p (x,q) — r{p,x) + r(p,q) — r(q,x), for any p, q, x in F. 

3. The discrete Laplacian L and the pseudo inverse L + 

To have a well-defined discrete Laplacian matrix L for a metrized graph T, we first choose 
a vertex set V(T) for T in such a way that there are no self- loops, and no multiple edges 
connecting any two vertices. This can be done for any graph by enlarging the vertex set by 
inserting more valence two vertices whenever needed. We will call such a vertex set V^(r) 
optimal. If distinct vertices p and q are the end points of an edge, we call them adjacent 
vertices. 

Note that a weighted graph corresponding to a metrized graph with optimal vertex set 
does not have self loops and multiple edges. Such graphs are usually called simple graphs in 
the literature. 

Let T be a graph with e edges and with an optimal vertex set V(T) containing n vertices. 
Fix an ordering of the vertices in V^(r). Let {L 1 ,L 2 ,--- , L e } be a labeling of the edge 
lengths. The n x n matrix A = (a pq ) given by 



ttpq 



if p = q, or p and q are not adjacent. 

if P 7^ q, and p and q are connected by an edge of length Lk 



is called the adjacency matrix of T. Let D = diag(d pp ) be the n x n diagonal matrix given 
by d pp = J2sev(r) a ps- Then L := D — A is defined to be the discrete Laplacian matrix of T. 
That is, L = {l pq ) where 



Ipq 



if p 7^ q, and p and q are not adjacent. 

if p ^ q, and p and q are connected by an edge of length Lk 
lps Hp = q 



l 



s6 v(D-{ p } 



The discrete Laplacian matrix is also known as the generalized (or the weighted) Laplacian 
matrix in the literature. 

Throughout this paper, all matrices will have entries in R. Given a matrix M, let M T , 
tr(M), be the transpose and trace of M, respectively. Let I n be the n x n identity matrix, 
and let O be the zero matrix (with the appropriate size if it is not specified). 
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A matrix M is called doubly centered, if both row and column sums are 0. That is, M is 
doubly centered iff MY = O and Y T M = O, where Y = [1, 1, ■ ■ • , f] r . 

Remark 3.1. For any graph T, the discrete Laplacian matrix L is symmetric and doubly 
centered. That is, l pg = 0, for each q G V(T), and l pq = l qp for each p, q G V'(r). 

P £V(r) 

In our case, T is connected by definition. Thus, L is of rank n — 1. Although L is not 
invertible, it has generalized inverses. In particular, it has the pseudo inverse L + , also known 
as the Moore-Penrose generalized inverse, which is uniquely determined by the following 
properties: 

z) LL+L = L, ii) L+LL+ = L + , Hi) (LL + ) T = LL + , iv) (L + L) T = L + L. 
The following properties hold for both L and L + : 

i) L and L + are symmetric, Hi) L and L + are EP matrices, 

ii) L and L + are doubly centered, iv) L and L + are positive semi-definite. 



For more information about L and L , see the article |C3] and references therein. 

I. Gutman and W. Xiao [I-W[ Lemma 3] obtained the following result when L arises from 
a graph where each edge length is 1 . 

Lemma 3.2. |D-Mt Equation 2.9] Let J be an n x n matrix having each entries 1, and 
let L be the discrete Laplacian of a graph (not necessarily with equal edge lengths). Then 
LL+ = L+L = I - ij. 



As an immediate consequence of Lemma 13. 2^ we have 
Corollary 3.3. Let V be a graph and let L be the corresponding discrete Laplacian matrix 
of size nxn. Then for any p, q G V(r), V" l+l sq = < " q ' 

The resistance distance has been studied in chemical literature and in computer science. 
See the articles |D-M] and |RBlj . The following fact plays an important role for the rest of 
the paper: 

Lemma 3.4. [RB2] [RBI J Let T be a graph with discrete Laplacian L and the resistance 
function r(x,y), and let H be a generalized inverse of L, i.e., LHL = L. Then r(p,q) = 
H P p — Rp q — Hg P + H q q, for any p, q G V(T). In particular, for the pseudo inverse L + we 
have r(p, q) = /+ - 2Z+ + /+ , for any p, q eV(T). 

Note that Lemma [3741 for the pseudo inverse L + follows from [D-Mt Theorem A]. Similarly, 
Babic et al [B KLNT] has Lemma [3.41 for L + and the graphs with all edges of length 1. 
Now, we can easily derive the following analogous result about the voltage function: 

Lemma 3.5. Let V be a graph with the discrete Laplacian L and the voltage function j x (y, z) . 
Then for any p, q, s in V(T), j p (q, s) = 1+, - Z+ - /+ + Z+. 

Proof. By equation (j2J), 2j p (q,s) = r(q,p) + r(s,p) — r(q,s). Then the result follows from 
Lemma 13.41 □ 
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4. Identities analogous to the extended Foster's identities 

Let r be a metrized graph with an optimal vertex set V(r) containing n vertices as before. 
Recall that if we consider T as a resistive electric circuit, the resistance of an edge is given 
by its length. If Lj is the length of an edge with end points p and q, then we have 
C pq = C qp = j- as the conductance of the edge e^. We write q ~ p when the vertex q is 
adjacent to the vertex p, i.e., when p and q are the end points of an edge. We set C pp = 
and C pq — if p and q are not adjacent, and define C p := C pq . It follows from the 

definitions that C p = l pp and C pq = —l pq if p ^ q, where L = (l pq ) is the discrete Laplacian 

of r. 

For any s G V(T), we |C1|. Section 5.5] proved the following identities by using Lemma [331 
and the properties of L and L + : 

P,gev(r) 
p~<?, p<q 

£ 3v,{q,sf WpCpq - V ~ 2 



C n 2 

w^jp^q, w<q 



p,q,WEV(T) V 



2 



E. / , \ C W pCpqCqt ^ — 3 1 C pq 



C n C n 2 2 CnCn 

p,q,t,weV(T) P q p,q€V(T) 1 ' 

w^pr^q^t, w<t P~9 

CnirtCmnC'inaCnt V 4 1 v \ ^nn 1 X ^ CriinCmnCn 



E. \ C up Cp W C W qC q t _ V — 4 1 \ (~y pq I 

3u{t,s) -— — + - 2^ T^ + o 2^ 



■"pw^wq^qp 



CnC in C n 2 2 ^ Cr.Cn 2 , , CrnC in C n 

p,q,t,u,w€V(r) p w q p,qeV{r) P q P) q,w&V{T) V w q 

u~p~w~q~t P~q p~g~to~p 

U<t 

Note that the formulas above are for graph paths containing at most 4 edges. A certain 
pattern involving conductances can be seen in these formulas. In the rest of the paper, we 
will extend formulas above, and obtain a new proof of the extended Foster's identities as a 
corollary. 

C 

Let P = (p it ) be the nxn matrix given by pu = — ^, and let P k = (p# ) be the fc-th power 



n C' c c 

of P. It follows from the definitions that p^ +1 ^ = %m l " lirri2 — — — whenever 

mi ,...,m k =l C i C ""i 1 ■ ■ • C mk 

n n 

k > 1. Moreover, we have T^Pa = 1 and that ^iPu = for any k > 1. As it is 

t=i i=i 
also given in the article [BCEGJ, P is known as transition probability matrix of the reversible 

Markov chain associated to T when T is considered as an electrical network, and P k is known 
as the k-step transition probability matrix. The value pj^ is the probability that the Markov 
chain attains vertex t at k-th step after starting from vertex i. 

n n 

Proposition 4.1. For any k > 1, we have £ Q/+p^ +1) = 1 - tr(P k ) + £ CJ+p^. 
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Proof. We first note that 

cvi ^7+<7- -/+/ — a i + - 1 n ^ <? ' 

i=l j=i 1 n 1 ~ 

The last equality follows from Corollary 13.31 For any k > 1 , we have 

E Cil itPu +l) = E °™ im2 E l tt C irnn Then by using equation 

i,t=l t,mi,...,m),=l i=l 



u i l itPit + 1^ c c n 2_, p -!-i 

;,t=l mi,...,m fe =l t=l mi rrtfc mi=l 

E + 1 E c ;r 2 ' ' '5*' - tr(P k ) 

i,t=l ^ i,mi,...,m*=l mi ' ' ' mfe 

n 1 n 

Y, C ^+- E Pl5t-tr(P k ). 



n 

i,t=l i,mi=l 



Then the result follows from the fact that p^ t = 1 = n for any fc > 1. □ 



t, mi=l mi=l 
n n 



We have E Cifet = E l t C « = " E + E C * 1 * = ~ n + 1 + E ^ where the 

i, f=l i,t=l i,t=l i i 

last equality follows from Corollary 13.31 Thus, the successive application of Proposition 14.11 
gives the following equality for any k > 1: 

n k—l n 

(4) E = k - n -H tr ( pl ) + E c *t- 

i,t=l i=l i=l 

Our main result is the following theorem which gives an identity analogues to the extended 
Foster's identities. Unlike extended Foster's identities which involve the resistance values, 
our formula involves the voltage values. 

n k—l 

Theorem 4.2. Let s e V(T). For any k > 1, we have E Ciji(s,t)p[^ = n — & + E tr(P'). 

i,t=l i=l 

Proof. Let s be a given vertex in V{T). For each k > 1, we have 

n n 

E C tJt (s, f)pg° = E C * ft* " # " E + &)P?\ b ^ LemmaEH 

i,t=l 

n n n 

- £) - E ^ - £)p^> since Epi fc) = L 



i,t=l i,t=l 



i=l i,i=l t=l 



= E ^ - E ^? ) > since E ^pl? = c - 

i=l i,t=l i=l 

Then the result follows from equation (jlj). □ 
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Next, we will state the extended Foster's identities as a corollary of Theorem 14.21 



^ n k—1 

Corollary 4.3. For any k > 1, we have - Cjr(z,t)p^ = n — k + tr(P'). 

i,t=l i=l 

Proof. By equation (JTj) , we have r{i,i) = ji(t,s) + j t (i,s) for any s 6 T. Thus, the result 
follows from Theorem 14.21 □ 

Note that the formula given in Corollary 14.31 was originally proved by E. Bendito, A. 
Carmona, A. M. Encinas and J. M. Gesto [BCEG} Proposition 2.3]. Next, we will use their 
methods to provide a second proof for Theorem 14.21 First, we recall some related definitions 
and results. 

Let e be the n x 1 vector whose entries equal to 1, and let e* denote the n x 1 z-th unit 
vector for each i = 1, . . . ,n. That is, = 1 if k = i and = if k ^ i. Suppose v % is the 
unique solution of the linear system Lu = e — ne l satisfying v\ = 0. The solution v % is called 
|BCEG] the equilibrium measure of the set V(T)\{i}. 

Since r(i,t) = -[v\ + v\) jBCEGl Proposition 2.1], we obtain the following proposition by 
using equation 

Proposition 4.4. For any i, j, t in ^(r), we have 

2 3l {s,t) = -(v l s + v s l+ vl + vl-v s t -^. 
On the other hand, we have [BCEGl Proof of Proposition 2.3] that 

. n k—1 

(5) -Y,C^=n-k + Y,^' 1 )- 

i,t=l i=l 

The second proof of Theorem 14.21 will be as follows. For any k > 1, we have 



2n 

i,t=i t,t=i 



n n 

CdiMv™ = — J2 + *f + A + »t- < ~ ^Vf, by Proposition 

i,t=l 

1 n n 

~ E C ^ > since C ^ = and E = 1 • 

i,t=l t=l 



k-1 

n — k + tr(P'), by equation (jSJ) 

i=i 
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